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1 Introduction 



During the last decade, starting with the pioneer paper by Minahan and Zarembo [1], a vast inte- 
grable structure has been unveiled in the N = 4 supersymmetric Yang-Mills (SYM) theory and in 
its dual supersymmetric string theory in the AdS§ x S 5 space-time. The integrable structure, known 
by the name of Integrability, was formulated in terms of an effective long range spin chain having 
PSL(2, 2|4) symmetry. For a recent review see [2]. 

The spectral problem, i.e. the problem of computing the anomalous dimensions of the gauge- 
invariant states in SYM and their string counterparts, is nowadays considered as conceptually solved. 
Spectacular progress was achieved in the last years in the computation of the gluon amplitudes [3,4] 
and Wilson loops [5, 6]. The next step towards the complete solution is to compute the three-point 
function of gauge-invariant operators representing traces of products of fundamental fields. This is 
obviously an extremely hard problem, but the encouraging developments over the last several years 
[7-11] raise the hope that Integrability can be used for this problem as well. 

Of special interest are the correlation functions of one-trace operators in the classical limit when 
the length of the traces is very large. Such operators are dual to extended classical strings in the 
AdS5xS 5 background. This correspondence allows to approach the problem both at strong and at 
weak coupling. The classical, or long-trace, operators are described in terms of algebraic curves, 
which appear as finite-gap solutions of the Bethe equations in the classical limit, or as classical solu- 
tions of the string sigma model [12-14] (see also the review [15]). 

On the string theory side, the problem of computing the correlation function of operators dual to 
classical spinning string solutions was addressed by several authors [16-22]. The three-point function 
of such operators can be thought of as a classical tunneling amplitude. To compute it, one should find 
the appropriate Euclidean classical solution for a world sheet embedded in AdS$ x S 5 and having 
the topology of a sphere with three punctures. However, there is not yet a consensus among the 
active workers on the field about the criteria to distinguish the relevant classical solution, neither 
there is an unambiguous prescription about how to construct the vertex operators associated with the 
punctures. The only case when the complete answer is known is that of two heavy and one light 
operators [17,23,24]. 

Alternatively, one can start with three one-trace operators in the weakly coupled gauge theory, 
where the computation of the correlation functions is a well defined problem. At tree level it is 
sufficient to count all possible planar sets of Wick contractions between the three operators. It was 
pointed out by Roiban and Volovich [8] that the calculation of correlation functions of gauge invariant 
operators reduces at the level of the spin chain to the calculation of the scalar products of states 
constructed out of B and C operators in the Algebraic Bethe Ansatz [25]. A systematic study of the 
case when the three operators belong to three different su(2) sectors of the theory was presented by 
Escobedo, Gromov, Sever and Vieira [9]. The tree-level correlation function of three su(2) operators 
was expressed in [9] in terms of scalar products of Bethe states in a periodic XXX spin chain with 
spin 1/2. The formalism developed in [9] was later applied to compute the correlation function of 
two heavy and one (or more) light operators [10,26], and the comparison with the strong coupling 
results on the string theory side [17, 23, 24] showed a precise match. The limit of three heavy, or 
classical, operators was then obtained in [11] in the case when one of the operators is BPS type, that 
is, protected by the supersymmetry. The main result of [1 1] is an elegant analytic formula in the form 
of a contour integral. The derivation was based on the expansion formula for the scalar products of 
two generic Bethe states due to Korepin [27]. 

In this paper we tackle the general case when neither of the three heavy operators is BPS type. We 
compute, at tree-level, the correlation function of three non-protected classical operators belonging to 
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three different su{2) sectors of the gauge theory. The principal object to be computed is the restricted 
scalar product of two Bethe states, in which part of the rapidities are frozen to a special value. Our 
result, which is a generalization of the main result of [11], was announced in the short note [28]. 
The analytic expression found in [28] is based on a factorization formula, which follows from the 
representation of the structure constant in terms of Slavnov-like determinants [29], proposed recently 
byFoda [30-32]. 

The plan of the paper is as follows. In Section 2 we review the basic notions about the XXX 
spin chain we are going to use, including the Gaudin norm, the Slavnov determinant formula for the 
scalar product, and its restricted version. Following Foda [30], we will perform the computations for 
the inhomogeneous spin chain, characterized by a set of external parameters (impurities) associated 
with the sites of the chain. Such a deformation of the problem allows to avoid ambiguous expressions 
containing poles and zeroes on the top of each other. In Section 3 we derive a representation of the 
Slavnov determinant in terms of free chiral fermions, and then perform bosonization. As a side result, 
we obtain a novel and potentially useful representation of the Izergin determinant for the domain 
wall partition function (DWPF) of the six- vertex model. In the limit when the number of Bethe roots 
tends to infinity, the bosonized expression decomposes into two computable factors, for which we find 
the analytic expression in Section 4. In section 5 we derive the structure constant of three classical 
non-BPS operators in the sl(2) sector of M = 4 SYM. 



2 Inner product of Bethe states in the inhomogeneous XXX chain 
2.1 The monodromy matrix 

The local fluctuation variable in a XXX spin 1/2 chain can be in two states, \, and f, which can be 
thought of as a basis of a two-dimensional linear space V on C. The spin chain is characterized by an 
isotropic Hamiltonian 

L 

H xxx = ( a m a m+l + a m a m+l + h a m a m+l) ■ C 2 - 1 ) 

m=l 

We will assume twisted periodic boundary conditions, 

<?t+L = K± a m > ( 2 - 2 ) 

which do not spoil the integrability, but allow to have a better control over the singularities. If the 
twist 

n = K-/K+ = e i * (2.3) 

is a pure phase, the twisted boundary conditions can be thought of as the effect of turning on a mag- 
netic flux of strength (f> [33, 34]. For us k will be an unrestricted complex parameter. 

In the framework of the Algebraic Bethe Ansatz [25], the spin chain is characterized by an R- 
matrix Rn{u, v) acting in the tensor product V\ <8> V 2 of two copies of the target space, 

R 12 (it, v) = u - v + iP 12 , (2.4) 

where P\ 2 is the permutation operator [35]. We will consider the inhomogeneous XXX spin chain, 
characterized by background parameters (impurities) 6 = {0\, . . . ,9l} associated with the L sites 
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of the chain. The twisted monodromy matrix T a (u) 6 End(V r ft ) is defined as the product of the 
i?-matrices along the spin chain and a twist matrix K = (q + «L) ' 

T a {u) = KR la (u, 6 1 + §i) R 2a (u, 9 2 + \i)... R La (u, 9 L + \%) = K . (2.5) 

For the homogeneous XXX spin chain all 6 m are equal to 0. The advantage of introducing the twist 
and the inhomogeneity parameters is that the expressions for some scalar products, which are ambigu- 
ous for 6 rn = and k = 1, becomes well defined for generic 6 m and k. 

The matrix elements A,B,C,V are operators in the Hilbert space V = V\ (g) . . . Vl of the spin 
chain. The commutation relations between the elements of the monodromy matrix are determined by 
the relation 

R\2{u — v)Ti(u)l2(v) = T2{v)T\{u) R\2{u — v) (2.6) 

which follows from the Yang-Baxter equation for R. As a consequence of (2.6), the operators B(u), 
C(u), and the transfer matrices 

Tin) = Ti a [T a {u)\ = A(u) + V{u), (2.7) 

form commuting families: 

[B(u),B(v)] = [C(u),C(v)} = [T(u), T(v)\ =0 for u, v £ C. (2.8) 

2.2 The Hilbert space as a Fock space for the pseudo-particles (the magnons) 

One can give the Hilbert space V a structure of a Fock space generated by the action of the operators 
B(u) on the pseudo-vacuum (t L | = | ft t )• The pseudo-vacuum is an eigenvector for the diagonal 
elements A and V and is annihilated by C: 

A(u) \t L ) = a{u) |t L ), V(u)\\ L ) = d(u)\t L ), C(u)\\ L ) = 0. (2.9) 

The dual Bethe states are generated by the action of the C-operators on the dual pseudo-vacuum 
(t L | = (tt 1 1> which has the properties 

(t L | A{u) = a(u) (| L |, (0|P(u) = d(u) (| L |, (t L \B(u) = 0. (2.10) 

The functions a(u) and d{u) depend on the representation of the algebra (2.6), while the i?-matrix 
(2.4) is universal. For the inhomogeneous XXX spin 1/2 magnet the functions a(u) and d(u) are 
given, according to (2.5), by 



a(u) = k + ]"[ (u - e m + i/2) , 

m=l 
L 

d{u) = k_ "[I ( u -°m- ■ 



m ~ l (2.11) 

L 



m=l 



The algebraic construction of the Hilbert space does not use the particular form of the functions 
a{u) and d(u), which can be considered as free functional parameters [27]. 1 The space of states is a 



'in this case one speaks of generalized SU(2) model. 



4 



closure of the linear span of all vectors of the form 



N 

|u> = n^K-)lt L ). (2-12) 

3=1 

Here and below we will use the shorthand notation 2 

u = u n = {ui,...u N }. (2.13) 

The operator B(u) can be viewed as a creation operator of a pseudo-particle (magnon) with momen- 
tum p = log • Similarly, the dual space of states is a closure of the linear span of all vectors of 
the form 3 

N 

<u| = (t L | Y[C( Uj ). (2.14) 

3=1 

Such states are generic, or off-shell, Bethe states. A Bethe state (2.12) which is also an eigenstate 
of the transfer matrix (2.7) is called on-shell state. For a chain of length L there are 2 L linearly 
independent on-shell states. Applying repeatedly the RTT relations (2.6) for the upper triangular 
elements of the monodromy matrix, namely 

A(v)B(u) = ^ B(u)A(v) - ^B(v)A(u), 
V{v)B{u) = 2=2=i B(u)A{v) + ^ v B(v)V(u), 

one obtains [35] the on-shell condition for the rapidities u = {ui, . . . un}, which coincide with the 
(twisted) Bethe equations 

g^TT "'--"* + * =-!, a = l,...,N. (2.16) 

Since the XXX Hamiltonian is hermitian, the set of roots u of a Bethe eigenstate must be invariant 
under complex conjugation. The corresponding eigenvalue T u (u) of the transfer matrix is 

Tu(u) = a{u) ^(ur + d(n) ^>T' (2 - 17) 

where Q u (u) is the Baxter polynomial 

JV 

Qu(u) = Y[(u- Uj ). (2.18) 

3=1 

The Bethe equations (2.16) can be considered as saddle-point equations for the Yang- Yang func- 
tional [25,36], which we denote by &~ u , 

8 Uj ^ = 2mnj (a = 1, . . . , N; rij £ Z). (2.19) 



2 We will indicate the cardinality only if this is required by the context. 
3 With the convention ^(u) 1 " = -C(u*) , the state dual to |v) is (-l) JV (v*|. 
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For the twisted periodic XXX^ spin chain, the Yang-Yang functional is given by 

N L 

%l = ^2 XI ( u i -°m + \i) ^g{uj -8 m + \i) - (uj -8 m - \i) \og{uj - 6 m - \i) 

j=\ m=l 

+ - u k + i) log(nj - u k + i) - (uj -u k -i) log(uj - u k - i)] (2.20) 

j<k 

N 

+ i<f> ^Uj. (2.21) 
3=1 

2.3 The pseudo-momentum 

For each set of points u = {uj}jL v define the function p u (u), called pseudomomentum 4 

e 2i Pu («) = Qu(u + i) = ^ Q (u-i/2) Q u (u + i) 
a(u) Q u (u — i) Qe(u + i/2) Q u (u - i) ' 

Here and below we will use the notation (2.18) for products over rapidities. In terms of the pseudo- 
momentum, the Bethe equations (2.16) read 



2iPuH = _i } u G u . (2.23) 



e 

The pseudo-momentum is determined by (2.22) modulo m. To characterize this function completely, 
it is necessary to specify a set of integers (mode numbers) {rij}^ 1 , not necessarily different, so that 

p{uj) = imj. 

2.4 Slavnov's determinant formula for the inner product 

The scalar product of two generic Bethe states can be computed from the algebra (2.6), the relations 
(2.9) and (2.10), and the convention {\ L \ \ L ) = 1. In general, the scalar product (v|u) of two 
Bethe states is given by a double sum over partitions of the sets u = and v = v^, which 
becomes increasingly difficult to tackle when number of magnons N becomes large. A significant 
simplification occurs when one of the two sets of rapidities, say u, satisfies the Bethe equations (2.16). 
It was discovered by N. Slavnov [29] that in this case the scalar product is a determinant. This is true 
for all integrable models with type R-matrix. 

Let the set u satisfy the Bethe equations (2.16). Then |u) is an eigenvector for the transfer matrix 
with eigenvalue T u (u), given by eq. (2.17). It was shown by Slavnov [29, 37] that the scalar product 
with a generic Bethe state (v| is proportional to the determinant of the matrix of the derivatives of 
T u (u), evaluated at the points of v, 

N 

(v|u) = Ha(v,)d( Uj )y u , v ; (2.24) 

3=1 

~ def 1 det jfc gf-T> fc ) 

n?=i a ( v j) det 3k ^r k 



4 There is no complete consensus about the terminology. The quantity we refer to as pseudo-momentum is related to the 

counting function Z{u) by e 2lpM = {-1) N+L e~ tZ( - u) . 
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(Eqs. (2.24)-(2.25) are equivalent to eq. (6.16) of [37].) The explicit expression for the matrix of the 
derivatives d Uj T u (v k ) is 

-dT u (v k )/duj = t( Uj - v k ) a{v k ) Q ^ Vk ~; } - t(v k - Uj ) d(v k ) + 3 

= a ( v k) ® U }; Vk ( tt(uj,v k ), (2.26) 

where the kernel Q,(u, v) is defined by 

n(u, v) = t(u -v)- e 2ip ^ v) t(v - u) , tin) = - — . (2.27) 

u u + i 

One can simplify the Slavnov determinantt (2.25) by noticing that ]lf = i = ^f^S^^i ■ 

The resulting expression is 

= det jk n( Uj ,v k ) 
detjfc u ._ Vk+i 

In order to get rid of the factors a(vj) d(uj) in (2.24), we rescale the annihilation/creation opera- 
tors as 

C(«)-XC(«) = ^, B(u)^M{u) = ^\ (2.29) 
a(u) a{u) 

and denote the rescaled inner product by (u|v): 

<viu> ^ (t L i n c (^) n B ( n )i t L ) = ^u,v . (2.30) 

Depending on the context, sometimes we will indicate the set of the inhomogeneity parameters and/or 
the cardinalities of the sets u, v and 6, 

(v|u> = <v|u>0 = (vjvlujvK. (2.31) 

2.5 The Gaudin norm 

The hermitian conjugation compatible with the scalar product (2.24) is C(u) = —B(uy, or 

C(u) = ^ B(«)t. (2.32) 
a(u) 

Assuming that the sets u and v are invariant under complex conjugation, u = u, v = v, and taking 
the limit v — > u in (2.28), one reproduces the determinant expression for the square of the norm of 
a Bethe eigenstate conjectured by Gaudin [38, 39] and proved by Korepin in [27]. The square of the 
norm is proportional to the Hessian of the Yang-Yang functional. In our normalization 

det jfc (igfe) „,„ 
(u|u) = - - 1 — . (2.33) 

The explicit expression for the matrix of the second derivatives is 

■ d*y u 2 ff 2 X- 1 ] 
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3 Operator factorization formula for the inner product 



3.1 Slavnov's determinant as a fermionic Fock space expectation value 

The Slavnov determinant expression (2.28) for the scalar product can be formulated in terms of free 
fermions and represents a tau-f unction of the KP/Toda hierarchy [40-42]. For the Izergin determinant 
this was shown and used in [31,43,44], see also the review paper [45]. The two sets of Toda times are 
related to the moments of the two sets of rapidities, u = {u±, . . . ,un} and v = {v±, . . . , vn}. 

The fermion representation we are going to use here is not the most natural one from the point 
of view of integrable hierarchies, but it reveals a hidden factorization property of the scalar product, 
which can be used to find an analytic expression in the thermodynamical limit N, L — > oo. 

Introduce a chiral Neveu-Schwarz fermion living in the rapidity complex plane and having mode 
expansion 

ij{ U )= Y, Au~ r ~K ij(u)= i>rU r ~l (3.1) 

rez+| rez+| 
The fermion modes are assumed to satisfy the anticommutation relations 

[VV, Vv]+ = [lp r , Vv]+ = , [ipr, VV']+ = 8 r ,r' , (3.2) 

and the left/right vacuum states are defined by 

(0|V>_ r = (0|Vv = and Vv |0> = $- r \0) = 0, for r > 0. (3.3) 

The operator -ip r creates a particle (or annihilates a hole) with mode number r and the operator ip r 
annihilates a particle (or creates a hole) with mode number r. The particles carry charge 1, while the 
holes carry charge —1. The charge zero vacuum states (3.3) are obtained by filling the Dirac see up to 
level zero. The left vacuum states with integer charge ±N are constructed as 

(N\ = l X (3.4) 

Any correlation function of the operators (3.1) is a determinant of two-point correlators 

(0|V>(u)^(«)|0) = ^-^. (3.5) 

Obviously, the Slavnov kernel (2.27) can be represented as the correlation function of two fermionic 
operators, located at the points u and v of the rapidity plane, 

Q(u, v) = (0| \i){v) - e 2ip i {v) i){v + i)\ [ip(u) - ij){u + i)\ |0). (3.6) 

The determinant of the matrix Q(uj,Vk) is equal to the correlation function of N pairs of such opera- 
tors, and the Slavnov inner product (2.28) can be given the following Fock space representation, 



<o| nf=i h Ktfj) - ^{uj + i)] uti gfa) - ^ {Vk) ^ + i)] |Q) 

<oinli^-+on£i^fc)io> 



^ni}=im + *niv=im;iu; 

(3.7) 
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Our aim is to rewrite (3.7) in a form convenient for taking the limit N — > oo. For that we first we 
transform the denominator, using the Cauchy identity, to 

N N 

£u,v = <oi n + *) n ^(«*)io> 
j=i k=i 

= det ' 



jk Uj - Vk + l 

A u A_ v 



N , (3-8) 



Ilj,fc=i(«j -Vk + i) 



Here and below we denote by A w the Vandermonde determinant associated with the set of complex 
numbers w = {w\, . . . , wn}, 



N 

def ■ 



A w = Y[(wj-w k ). (3.9) 

j<k 

Then we represent the expectation value in the numerator of (3.7) as a product of difference operators 
acting on the Cauchy product (3.8). The result is 

1 N 

y uv = _ Yl (l - e 2ip «to> e i9/9v ^ (e- id / 8u i - l) /C U)V , (3.10) 
u,v j=i 

where e l9 ^ 9u denotes the shift operator u — > u + i. We can commute the denominator of the Cauchy 
product to the left, using the relations 



N N 

N N 
e id/dvj TT = E-( Vj )T\ jd/dvi, (3.11) 

until it cancels its inverse. Here and below we denote 

± Q u (v±i) ± Q v (u±i) 

^ (v) = ^r' ^ (u) = -^r (3 - 12) 

Now we can write (3.10) in a factorized operator form, 

N 



v i=i 

r J] (l - E+( Uj ) e- id ' du i) A u • 1. (3.13) 



The factorization is not complete because we must commute the operators e ld / dv ^ to the right through 
the functions in the second factor, which depend implicitly of the v-variables: 

e"/** £+[«,-] = {l- { l k)2 + 1 ) E t\uj] (3.14) 
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3.2 Another writing of the operator factorization 



The operator representation (3.13) of the inner product of an on-shell state ||u ) and an off-shell state 
(v|| is the main tool we are going to use to investigate the classical limit of large L and N. Here 
we will give it a more abstract formulation in terms of a pair of non-commuting functional variables, 
which become c-functions in the classical limit. 

For any set of points u = {uj}^ =l in the complex plane and for any function f(u) define the 
functional 

1 N 

def 1 



<[/] = ^n(w(^) e±<8/8wi ) Au ' (3.i5) 

or, in terms of free fermions, 

, ±m (N\Uk=i[^k)-f(u k )^(u k ±i)]\0) 



/ \f] 

wn£=i#«k)io> 



det jfc (u*" 1 -/(«,-) K-±0 fc_1 ) 



det 



jk 



(-r 1 ) 



(3.16) 



The functional srf^ [/] are completely symmetric polynomials of the values of the function / on the 
set u, having total degree N. 

Assuming that u n v = 0, we can write the r.h.s. of (3.13) as a matrix element 

_ » (j^£g , (3 , 7) 
(v|u) 

where the operator functions U (v) and V(u) satisfy the algebra 

U(v)V(u) = V(v)U(u) (l - (M _J )2 + 1 ) (3-18) 

and act on the vectors (v| and |u) as 

U(v) |u) = e 2ip ^ E~(v) |u) = ^E+(v) |u) , 

a(v) u (3.19) 

(v\V(u) = E+(u) (v|. 



In this way the problem of evaluating the inner product reduces to the problem of evaluating the 
functionals (3.15). 



3.3 Generalization to non-highest-weight states 

The A-magnon states (2.12) are highest weight states with respect to the fully ordered state, pseudo- 
vacuum | | L ). Each operator M(uj) flips one spin down so that the third component of the spin of the 
state with N magnons built on the pseudo vacuum |0) is = — A^. 

A complete system of states is obtained by acting with magnon creation operators on the vacuum 
descendant states (§~)^| | L ) with K < L/2, which correspond to ferromagnetic vacua rotated away 
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from the third axis. On the other hand, the components S 1 * 1 of the total spin can be obtained by taking 
the infinite rapidity limit of the magnon creation and annihilation operators, 

B(u)~-S - , C(u)~-S+. (3.20) 
u u 

The factor l/u is obtained by comparing the large u asymptotics W(u)M(u) | t L ) — L/u 2 , which 
follows from (2.34), with the normalization of the spin operator (t L |S + S~| | L ) = L. Therefore an 
M-magnon state built upon a vacuum descendent can be obtained as the limit of a iV-magnon state 
(2.12) with K = N - M of the rapidities sent to infinity: 

N-K' 

<t L i(s + ) x ' n %)i=(v M ,u^'i, 

j=1 (3.21) 

N-K 

H B(itjfc) (S-)^|t L ) = |u A r_ i ,Uoo^), 
1=1 

where the limits should be taken sequentially according to the definition 

Kr-iUoo)^ lim ^|u JV )=S-|u JV _i>. (3.22) 

u N — >co i 

To evaluate the inner products of such states, we need to compute the result of sending to 
infinity in the functional ^^[f^], assuming that the function / ± (u) behaves at infinity as 

f^u) ~ e TiK±/u , u^oo. (3.23) 

From the definition (3.15) we find, taking into account that A UiV ~ u^ r ~ 1 A UN _ l , 



= lim _ e TiK/u Ne ±id/du N)u N-l A ^ 

MJV^OO A UJV _ 1 

= ±(K ± -N + l)At N _ 1 [f ± ]. (3.24) 

Applying this relation the necessary number of times, one can obtain the generalization of the operator 
factorization formula (3.17) to the case of an inner product of non-highest-weight states. 

3.4 The case when the Bethe eigenstate is a vacuum descendent 

The inner product of Bethe state with a vacuum descendent is obtained from eq. (3.17) by sending 
all the roots from the set u to infinity. In this limit limit the sets u and v are infinitely separated, 
the commutator in (3.18) vanishes, and (3.19) gives U(v) = and V(u) ~ e *JV/« Applying 
sequentially (3.24) to all variables u, one obtains srf~ N [V] = (-1) N N\ , and (3.17) gives 

N 

(t L | n C ^') i^) N \t L ) = ^oc- VJV = N\ K + [d/a}. (3.25) 

To evaluate the scalar product of two vacuum descendants, we have to send the set v to infinity as 
well. Applying sequentially (3.24), with if_ = L, to all variables from the set v, one finds 

(t i |(S + ) iV (S-) JV |t L ) = ^{oc-Moc-} = (Nl) 2 • (3-26) 
The second factor counts the number of ways to have N reversed spins in a chain of length L. 
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3.5 Two-kink pseudo- vacua and restricted scalar products 

In view of the applications to S YM, we are going to consider pseudo-vacua with two kinks at distance 
K. In such a state the first K spins are down and the rest L — K spins are oriented up: 

(i K t L ~ K \ = (U ■» 4-TT - L ~ K 1 1 • (3-27) 

This state can be obtained by acting on the left pseudo- vacuum (| L | with K annihilation operators 
[30] 

K 

{\- k \ l - k \= <t L in c (*i)=< z *i' ( 3 - 28 > 

with rapidities z, K = {z\, . . . , z k } determined by the values of the inhomogeneity parameters associ- 
ated with the first K sites: 

z k = 6 k + \i (k = l,...,K). (3.29) 

Hence, the inner products with the left pseudo- vacuum replaced by the two-kink state (3.27) are 
evaluated by restricting K of the rapidities 5 

N—K N 

(i K t L - K \ [J C( Vj )l[M(u k )\t L } = (v^UzJu^, (3.30) 

j=l k=l 

with 

e L =e L _ K ue K , x K = e K -\i. (3.3i) 

The restricted scalar product (3.30) can be computed by a factorized expression similar to (3.17), 

(vuzlu). = Hf ^MBH , (3.32, 

where the operator functions U(v) and V{u) satisfy the algebra (3.18) and act on the vectors (v| and 
|u) as 

U(v) |u) = K 1 E^v)\u) , 

EJ (v) Qe(y + p) (3.33) 

(v\V(u) =E+(u)E+(u) (v|. 

This trivial substitution shows the power of the operator expression (3.13). As a comparison, when 
evaluating the restricted scalar product using the original Slavnov determinant, one comes upon spuri- 
ous singularities (poles and zeroes on the top of each other), which require multiple use of l'Hopital's 
rule. 

In the limit u, v — > oo, defined as in (3.24), one obtains 

^oc- oc;v-k Uzk = N\(N - K)\ ( L ~ . (3.34) 

The second factor in (3.34) counts the number of non-equivalent ways to reverse N spins among the 
remaining L — K up-spins of the partially ordered pseudo-vacuum. 



5 The restricted inner product (3.30) has been studied in [30,31,48,49]. A statistical interpretation of the restricted scalar 
product as partition functions of the six- vertex model is given in [31,49]. 
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3.6 Gaudin-Izergin determinant and pDWPF 

In the particular case K = L = N, the restricted scalar product 

N 

fu,z = k|u> w =(fll]%) It"), z N =9 N + i/2, (3.35) 

3=1 

evaluates the partition function of the six- vertex model with domain-wall boundary conditions (DWBC) 
on a iV x N square grid [27,50]. With this specialization of the rapidities, the second term of the kernel 
£l(u, v), eq. (2.27), vanishes at v = Zk and the Slavnov formula (2.28) gives 

The determinant representation (3.36) of the DWBC partition function was obtained by Izergin [51, 
52]. For the first time the ratio of determinants (3.36) appeared in the works of M. Gaudin [38,53] as 
the scalar product of two Bethe wave functions for a Bose gas with point-like interaction on an infinite 
line. 

The Gaudin-Izergin determinant (3.36) has an obvious writing in terms of free fermions, 

x _ jgl nf = i Muj) - M*, + Q] nti #**)|0) 33 

(oinf = i^K+onf=ife)io) 

and is expressed in terms of the functionals (3.15) as 

2f u , z = (-l) N a?-[E+]. (3.38) 

The limit of iF UjZ when part of the rapidities u are sent to infinity was recently studied by Foda and 
Wheeler in [47] and given the name partial domain-wall partition function, or pDWPF. We compute 
the pDWPF by applying sequentially (3.24) to N — n of the variables, with the result 

^u„uoo--n iZjv = (-1)" (N - n)\ zf-jE+J . (3.39) 

An alternative proof of eq. (3.39) is presented in [47]. The fermionic representation of pDWPF is 



v , (N - n\ U]=i Muj) ~ ^ + »)] Uti 1>{zk)\0) 
uoo N -™ z*r = (N-ny. T7 — = . (3.40) 

" <^-nin?=i^(«i+on^i^*)io> 



The expectation value on the r.h.s. of (3.40) is defined for any pair of non-negative integers N and n, 
but it vanishes identically when n > N, This yields a pair of identities 



< ± J^z F J=0 for N<n. (3.41) 

3.7 Properties of the functionals [/] 
Expansions 

The functionals s^^, defined by eq. (3.15), are obviously completely symmetric polynomials of 
degree iV of the variables f(ui), . . . , /(itjv). The coefficients of the polynomial are obtained by 
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expanding the product in (3.15) as a sum of monomials labeled by all possible partitions of the set u 
into two disjoint subsets u' and u", with u' U u" = u, 

< = e f n /<«')) lr n ^ e/au 'A u . 0.42) 

u'Uu"=u Vu'eu' / U u'Su' 

Here |u'| stands for the number of elements of the subset u'. The last factor is evaluated as 

1 / n^ a/8 "'Vu= n "^i^P- 0.43, 



\u'&u' / u'Gu'X'eu" 



The expansion (3.42) can be used as an alternative definition of the functional =e^[/]. For f = d/a, 
this expansion was thoroughly studied by Gromov, Sever and Vieira [1 1]. It was found in [1 1] that for 
constant function f(u) = n, the expansion (3.42) does not depend on the positions of the rapidities u 
and the functional srf^ [f] is given in this case by 



= (1 - k) n = exp ( \ J ■ (3-44) 



n=l 



2. The linear term in / as a contour integral 

The linear term in / can be evaluated as a contour integral: 

N 



j^f{u)Et{u) + 0[f]. (3.45) 



= 1 ± 

The integration contour A u encircles all points of the set u and leaves outside the possible singularities 
of the function f(u). 

Functional identities 

Using the fermionic representation (3.16) and the fact that the fermion correlator is translation invari- 
ant, we transform 

(N\ Uk=i $(vk-i)-f(vk)$(vk)]\Q) 
(N\Uk=i^k)\0) 

(0\ULi\^k) - j^^v k -i)]\N) 

= (-l) N f( Vl )...f(v N ) 



(0\UtMvk)\N) 

= (-1) N /(«!)... /KX"[1//]. (3.46) 
Hence, £? + and are related by the functional identities 



<*[!//] = (-1)^ =^ ~ [ !} ; • (3.47) 

IL=1 /(«; 
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Figure 1: An example of the distributions Figure 2: The cuts C u and C v and the inte- 

u and v for N = 50, each consisting of a gration contours A u and A v for the one-cut 

single macroscopic Bethe string. solution of Fig. 1. 



4 Classical limit 

In this section we will find the classical limit of the inner product (2.28). The classical limit is achieved 
when L, N — > oo with a = N/L finite, and some additional assumptions on the distribution of the 
rapidities ujv and vjy. In the condensed matter literature the classical limit, in which each Bethe 
string has macroscopic number of particles, has been studied by Sutherland [34] and by Dhar and 
Shastry [54]. In this regime the roots u condense along a curve C u in the rapidity plane, consisting of 
several connected components C Ul , . . . , C Un , with ui U . . . u n = u, symmetric about the real axis, 
with slowly varying linear density p u (u) [12]. The curve contains Nf. = particles, 

/ p{u)du = N k , N l + --- + N n = N. (4.1) 

We assume that the filling fractions = N^/L associated with the cuts C k remain finite when 
L — > oo. Then the size of each curve is ~ L. We assume a similar behavior for the rapidities v. An 
example of distributions u and v with = 50 and n = 1 is given in Fig. 1 . 

In the classical limit, the arguments U and V in Eq. (3.17) become c-functions, and the inner 
product factorizes to 

y uw = (-l) w <+[ Ke iG "- !G «]<[e ,G '], (4.2) 

where 

G u (u) = d u logQ u (u), G v {u) = d u \ogQ v (u), G e (u) = d u \ogQ s (u) (4.3) 

are the resolvents associated respectively with the sets u, v and 6. The resolvent G u (u) is a meromor- 
phic function of u with cuts C„, . . . , C™ and asymptotics N/u at infinity. The discontinuity across 
the cuts is proportional to the density p u {u). 
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The form of the curve C u and the density of the distribution of the roots is determined by the 
finite-zone solution constructed in [12], which we discrete shortly below. In the vicinity of each cut 
C„, the pseudo-momentum 

p(u) = G u (u) - \G e {u) + \<\> (mod tt) (4.4) 

splits into a continuous part i>{u), equal to the half-sum of the values of p{u) on both sides of C u and 
a discontinuous part p(u), which is proportional to the density p u : 

Pu(«) =/u(«)+Pu(«), \Pu(u)\ = 7rp u (u). (4.5) 

When the set u satisfies the Bethe equations, the eigenvalue of the transfer matrix T u = 2 cosp u is 
analytic in u and hence takes the same value on both sides of the cuts C„. This yields the boundary 
condition sinj5 u sin j$ = on the cuts. Along each cut p u > 0, which implies ^ u = (mod tt), or 

= 7m k , u e C* (4.6) 

where n k is the mode number associated with the cut C„. The branch points are zeroes of the entire 
function A(u) = T^(u) — 4 = — 4sin 2 p u . The forbidden zones A(u) > are associated with the 
cuts of the pseudo-momentum p(u) on the first sheet. 

The typical situation in the homogeneous limit is when A(u) has a double zero at u = oo, 2n 
simple zeroes a±, a±, . . . , a n , a n , and infinitely many negative double zeros at a_i, a_2, • • . a_fc, . . . 
where p(a_^) = 27rn_ / t, accumulating at the point of essential singularity u = 0. The cuts are along 
the forbidden zones between and (k = 1, . . . , n). 

The derivative of pseudomomentum dp u (u) as well as the exponential e 2ip<M ) are defined on a 
Riemann surface associated with the hyper-elliptic complex curve 

n 

V 2 = \{{n 2 -al). (4.7) 

k=i 

The values of e tp ^ on the first and on the second sheet are related by 

The derivatives of the pseudo-momentum in Nj, 

d Nj p(u)du = ujj(u), j = l,...,n (4.9) 
form a basis of holomorphic abelian differentials on the complex curve curve: 

^— f ujj = 5 kj , k,j = l,...,n. (4.10) 
2m J A k 

The cycle A 1 ^ represents a closed contour on the Riemann surface which encircles the cut anti- 
clockwise, as shown in Fig. 2. 

By the factorization formula (4.2), the computation of the inner product in the quasiclassical limit 
reduces to that of the functionals s^^. This last problem was solved in the particular case f = d/a 
in [11]. Below we give an alternative, and a little bit shorter, derivation of their result. Let us remark 
that although both derivations seem rather convincing, neither of them is a rigorous one. 
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4.1 Classical limit of the functionals srf^- [/] 

Assume that the rapidities u are distributed along a (possibly milti-component) contour C u in the 
complex u-plane with finite and sufficiently smooth density When the functional argument / is 
small, by eq. (3.45), 



io g <m = ±yg e ^w + ou% 



(4.11) 



where the integration contour ^4 U encircles C u anticlockwise, and the function q ± (u) is defined by 



q ± (u) = -i\og[f(u)]±G u (u). 



(4.12) 



On the other hand, the functional relations (3.47) allow us to determine the asymptotics for large /, 
which can be written again in the form of a contour integral. To see that we first express 



log 



N 



(-i) N Hf( Uj ) 

3=1 



= <f^,G u (u)(log\f(u)]+m) 



(4.13) 



Substituting (4.13) in (3.47), we find the large / asymptotics 

log^l/l^^i^M+wfTe-^W) + 0(/- 2 ]. 



(4-14) 



We will look for a solution compatible with the behavior at / — > 0, oo, which should be of the form 

\og^[f}=l ^(e^M), (4.15) 
where q(z) is defined by (4.12) and the meromorphic function F(co) has asymptotics 



' ±uj + 0{oj 2 ) if \u\ <C 1, 

Tl log(-w) 2 T 1/w + 0{l/u 2 ) if |w| > 1. 



(4-16) 



The function F(u) can be determined completely by comparing the Ansatz (4. 15) with the known 
exact solution (3.44) for / constant. Assume that the function F ± (uj) is expanded in a Taylor series 



F±(u) = Y,Fn" n 



(4.17) 



n>l 



in some vicinity of the point to = and compute the r.h.s. of (4.15) for q ± (u) = —i log k ± G u (u), 
which corresponds to f(u) = k. The contour integral can be evaluated by expanding the contour to 
infinity, and we find 



£ F - f t e±mq[u) = e F - f t Kn (i ± 4 + • • • )n = t s f « < 4 - 



18) 
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Comparing (4.18) and (3.44) we find that F± = ±l/n 2 and that the Taylor expansion (4.17) is that 
of the dilogarithm, 



id 

F ± (oo) = ± -2 = ± Li 2(w). (4.19) 

n=l ^ 

The asymptotic behavior (4.16) is satisfied thanks to the functional equation for the dilogarithm, 

Li 2 (-) = -U 2 (w) ~^r-\ log 2 (-u;). (4.20) 

Moreover, the property (4.20) of the dilogarithm leads to a pair of functional equations for £^[f], 
which are the scaling limit of (3.47). 

If the resolvent G u has several cuts on C„, . . . , C™, then the integration contour in (4.15) splits 
into n disjoint contours A^, . . . , A™, and the functional srf^ [f] is given in the classical limit by 

log*C[/] ~ ±y^U 2 (/(«)e ±iG »M), A u = Ug =1 A*. (4.21) 

The fc-th term grows as a^L, where ctfc = N^/L is the filling fraction associated with the cut C„. Let 
us emphasize that in the derivation of (4.21) we did not assume that the set u is on-shell. 

Finally, let us note that the functional identity (3.47), or equivalently, the property of the diloga- 
rithm (4.20), leads to a second integral representation, 

log< ± [/]^T^Li 2 (/- 1 (n) e ^ G ^))+£ ^G u («)log/(«) +Z7TN. (422) 



4.2 Classical limit of the Slavnov inner product y u , v 

Substituting (4.22) in (4.2), we write the logarithm of the scalar product as a contour integral 

log ,5V = zvriV + j> ^ Li 2 ( Ke iG u (u) + iGAu)-iGe(u)^ _ j dM ^ ^_ iGu(tt)+iGv(tt) ^ (423) 

The r.h.s. of (4.23) can be reformulated in terms of a contour integral around the ensemble of the cuts 
of the function 

q(u) = G u (u) + G v (u) - G e {u) + (4.24) 

on the physical sheet of its Riemann surface, with integrand depending only on q(u). This follows 
from the fact that the resolvent G u satisfies on its cuts the boundary condition (4.6), 

2G u (u) -G {u) + d> = 2im k for u £ C„. (4.25) 

(Here G u is the half-sum of the values of the resolvent on both sides of and n k € Z is the 
corresponding mode number.) Hence, if is the value of the function q(z) on the physical sheet, 
defined by (4.24), then the value of q{u) on the second sheet under the cut is given by q^ = 
—G u + G v . We conclude that the two integrals in (4.23) have the same integrand Li 2 (e* 9 ), but the 
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Figure 3: The contour A v and the 
deformed contour A^ for the integral in 
(4.26). The contour A™ starts at z = — ioo 
on the second sheet, passes on the first sheet 
at the lower branch point, continues to the 
upper branch point, where it returns to the 
second sheet and continues to z = +ioo. 



Figure 4: The contour A u for the integral 
in (4.27), obtained as the limit v — > u of the 
integral (4.26). When C v -> C u , the two 
logarithmic branch points on the first sheet 
join the two simple branch points at the ends 
of the cut C u . 



contours A u are placed on the second sheet of the Riemann surface of the function q(u). After pulling 
all connected components A 1 ^ C A u up to the first sheet across the cuts C*, eq. (4.23) takes the form 

log^ u ,v = wrJV+ j> ^Ii 2 (e iff(u) ). (4.26) 
A U UA V 

(The minus sign is compensated by the change of the orientation of the contours A^ after they are 
moved to the first sheet.) 

The integral along A^ is however ambiguous, because the integrand has two logarithmic cuts 
which start at u = oo on the second sheet, cross the cut and end at two branch points on the first 
sheet. The ambiguity is resolved by deforming the contour A u to a contour A^ 3 which encircles also 
the point z = oo on the second sheet. 6 In the case of a one-cut solution, the contour A™ is depicted 
in Fig. 3. 

In the limit u — > oo, eq. (4.26) must reproduce the result of [11] about the scalar product of two 
general Bethe states and a vacuum descendent (eqs. (3.28)-(3.29) of [11]). This is indeed the case. In 
the limit u — > oo, the integration goes only along the contour C v and the function q in the integrand 
is given in the homogeneous limit by q = G v — ~ + (f>. 

4.3 Classical limit of the Gaudin norm 

An expression for the square of the Gaudin norm in the classical limit can be formally obtained from 
(4.26) by taking G v = G u . Here we assumed that the two sets of rapidities are invariant under 
complex conjugation, u = u and v = v. When v — > u, the integration contour A u in (4.26) can be 

6 The author is indebted to Nikolay Gromov for performing a numerical test and for suggesting how to place the integra- 
tion contours. 
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closed around C u and C v as in Fig. 4, and in the integrand one can replace q(z) — > 2p(z), where p(u) 
is the quasi-momentum (4.4). We find for the square of the Gaudin norm (u||u) 2 = =5^ U)U , with 

log^u,u = f^U 2 (f*M). (4.27) 

The term iirN is compensated by another such term which appears because of the hermitian conju- 
gation ,6^ = —C. Furthermore, when = C^, the two cuts on the second sheet end at the two 
branch points where 1 — e 2 * Pu = 2ie LPu sinp u = and Li2(l — e 2iPu ) has a logarithmic singularity 
L12 ~ sin(p u ) logsin(p u ). Therefore there is no obstructions for placing the integration contour A u . 

An expression of the Gaudin norm as a linear integral was derived in [11]. One can check, using 
the fact that p u (z) = ±iirp u (z) on the two edges of the cut, that the contour integral (4.27) can be 
transformed into (twice) the linear integral in eq. (2.15) of [11]. 



4.4 Classical limit of the restricted inner product 

From Eqs. (3.32)-(3.33), one readily evaluates the classical limit of the restricted scalar product (3.30), 
y u ,vuz = (-1) N K + u z lKe iG »- lG °- iG *] <-[e^ +iC H (4.28) 
where U z are the inhomogeneity parameters, and then use (4.21). We find 

log^ u ,vUz = mN+j^ Li 2 (^e iG ^ +iG ^- iG °W)-j ^ Li 2 ( e -iGu(«)-HGv(«)-HG.(«)) 
Av An 

= i7T N+ j> ^u 2 (e lG ^ +lG ^- lG °^y (4.29) 
A U UA V 

In the second line we used the classical Bethe equation for 4> = 

2f u (u) - G e {u) - G x {u) = 0(mod2vr), u G C u . (4.30) 

We see that the restricted scalar product ^ UiV uz is given by the same contour integral (4.26), where 
the integrand depends only on u, v and 6. The dependence on z is through the boundary condition on 
the cuts of the resolvent G u . 



4.5 The derivatives in N k 

Knowing the hyper-elliptic Riemann surface (4.7), it is possible to compute the logarithmic derivatives 
of the scalar product with respect to the filling numbers N^. For that we use the fact that, according to 
(4.9), the derivatives of the pseudo-momentum form a basis of abelian differentials for the Riemann 
surface. 

Take for example the Gaudin norm. From (4.27) we find for its logarithmic derivative 

<91ogXi,u . / dz dp u 
= —Zi 



fdz&pv / MpuW N 

- u 

^ j> u k (z) log(2sinp u (z)) . 



(4.31) 



The differentials 0Jk{u) = dN k Pu{u)du are analytic on the hyper-elliptic curve (4.7) and are obtained 
as a solution of the relations (4.10). 
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5 3-point functions of trace operators in Af = 4 SYM 



5.1 The three-point function 

In a su(2) sector of the SYM theory, the operators are made of two complex scalars. We consider, as 
in [9-11], the correlation function of three single-trace operators of the type 

d ~ tt[z Li - Ni x Ni + ...], 

2 ~Tt{Z L2 - N2 X N2 + ...], (5.1) 
3 ~Tt{Z L3 - n *X Ns + ...], 

where the omitted terms are weighted products of the same constituents taken in different order. The 
weights are chosen so that the operator Oi is an eigenstates of the dilatation operator with dimensions 
Aj. The two-point and the three -point functions are determined, up to multiplicative factors, by the 
conformal invariance of the theory, 

( OiixJOjixj) > = U Stj Nl (5.2) 



/ n i \n ( \n i \\ L ± L 2 L 3 y 'NiN 2 N 3 Ci 23 (A) 

( 1 (x 1 )0 2 (x 2 )0 3 (x 3 ) > = |xi2|Al+A2 _ A3|x23|A2+A3 „ Al|:r3i|A3+Al _ A2 , (5-3) 

where Mi are arbitrary normalization factors. 7 The factor Lj count for the cyclic rotations of the trace 
operator Oi. The structure constant Ci23(A) has perturbative expansion 

N c C 123 (X) = c{% + XCill + ... , (5.4) 
where N c is the number of colors and A is the 't Hooft coupling. 



5.2 The structure constant in terms of scalar products of Bethe states 

At tree level, the structure coefficient is a sum over all possible ways to perform the Wick contractions 
between the fundamental fields and their conjugates. A non-zero result is obtained only if 

iVi = N 2 + iV 3 (5.5) 

and the number of contractions Ly = \(Li — Lj — Lk) between the operators Oi and Oj are 

L 12 = Li - N 3 , L 13 = N 3 , L 23 = L 3 - N 3 . (5.6) 

The product of all free propagators in the contractions between Oi and Oj reproduce the factor 
\xij\~ Ai ~ Aj+Ak in (5.3), with tree-level conformal dimensions A, = A^ = Li. By planarity, 
all Z fields in 3 must contract with Z fields in 2 and all X fields in 3 must contract with X fields 
in 0i, as shown in Fig. 5. Hence there is a single term in the sum in Eq. (5.1), Tr(Z L ™X Ns ), for 
which the planar contractions with 0\ and 3 do not vanish. 

In order to compute the tree-level structure constant cf 2 \ by the method of [9], one should know 
the wave functions at one loop. At one loop level, the operator Oi is represented by a /Vj-magnon 

7 Our normalization is slightly different that the normalization used in [9], namely A^ erc = A/? crc / 'Li and Cjf| = 

Cliffy V £1X2X3. 
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Figure 5: An example of a planar set of contractions contributing to the three -point 
function of the operators 0\, O2 and O3. The Z — Z and the X — X propagators are 
represented respectively by continuous and dashed lines. 



Bethe eigenstate 8 with energy Aj of the periodic XXXi/ 2 spin chain of length Lj (i = 1, 2, 3). Let us 
choose the pseudo vacua in the three chains as 

Tr(Z Ll ) -> I | Ll ) , Tr(Z L2 ) -> || f L2 ) , Tr(Z L3 ) -> | | L:i ) . (5.7) 

Then the three operators are determined by three Bethe eigenstates 

0i -Hu)^, 02^||v) i2 , 3 ^|w) L3 , (5.8) 

with u = {tii, • • • ,unA, v = {ui, . . . , v at 2 }, w = {u>i, . . . , u>at 3 }- A natural normalization of the 
two-point function is given by choosing 

M = (u|u>, AT 2 = (v|v), AA 3 = (w|w). (5.9) 

The tailoring (= cutting+flipping+sewing) prescription of [9] gives 



,10) = (v||u) Ll2 |w) Ll 

V'( U l U )ii< V ll V )i2H w )L3 



CI23 = 7==== == ■ (5-10) 



Since the inner products in the numerator are evaluated for subchains, the Bethe vectors are no 
more on shell. However, as proposed in Ref. [30], one can replace (v||u) l 12 by an inner product of the 
type (. . . lu)^, which can be evaluated by the Slavnov formula. This can be done after deforming the 
problem by introducing impurities with rapidities = at the sites of the i-th spin chain 



8 For simplicity consider highest-weight states. The generalization to arbitrary states is outlined in section 3.5. 
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(i = 1, 2, 3). We denote the rapidities associated with the contractions between the operators Oi and 
Oj by 0fa), so that = #( 12 ) U 0( 13 \ etc. Along the lines of [30] we obtain for the structure 
constant, up to a phase factor, 

c (0) = JvUzJu^q) (z|w) g( 3) ^ z = 0{13) + . /2 _ (5 n) 

V < u l u )e(D ( v l v )e(2) < w ll w >0(3) 

In this way we expressed the structure constant in terms of the quantities evaluated in Sections 3.5 and 
3.6. 

5.3 The BPS limit 

The structure coefficient should be normalized so that in the limit when all rapidities go to infinity it 
tends to the structure coefficient for three BPS fields [55], 

lim C[% = Cff. (5.12) 

u,v,w— >co 

From (3.34) we obtain the expected result 

123 = / = ?= = = ■ (5.13) 



f'' , ft) ,w (».)" w 



,2 



L 3 \ //LA /LA /L 3 

^3/ V \nJ \nJ \N 3 



5.4 The limit of three classical operators 

Here we take the limit when the three fields become classical, 

Li^oo with ai = ^- fixed (i = 1,2,3). (5.14) 
Lj 

using the results of Section 4. The two factors in the numerator in (5.11) are evaluated using (4.29), 

/du 
^Li 2 (expi[G u + G v -G 0(12 )]) (5.15) 

A U UA V 

/I 
— Li 2 (expi[p u + p v + lG ei3) ]) , 

A U UA V 

/du 
— Li 2 (expi[G w - G 0( 23)]) (5.16) 

A w 

/du 
— Li 2 (expi[p w + \G B{ 2) - \G eW }) , 

where we introduced the three pseudo-momenta (for = 0) 

Pu = G u - \G 6 (i), p v = G v - ^G ( 2 ), p w = G w — \G B (z). (5.17) 
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The norms in the denominator are evaluated using (4.27). Collecting all terms we find 



- I / t ^ (^") " I / U * ( e " Pv ) " § J t U * ( e2 " W ) • t 5 " 18 ) 

Aw 

The tree-level structure constant is obtained by setting all inhomogeneity parameters to zero: 

log C$3 ~ j> ^U 2 {e ip ^ +l ^ +lL ^ 2u ) + £^u 2 (e* p ^ +l ^- L ^ 2u ^ 

A U UA V A m (5.19) 

-\f p Li 2 (e 2 ^) - \ f p Li 2 (e 2 ^(«)) - § /* ^ Li 2 ( e *M«)) . 

To make connection with the result of [11], one should send all u's to infinity, which is the same as 
taking C u = and neglecting the integration along ^4 U . 

6 Conclusions and speculations 

The principal result of this work is the operator factorization formula (3.17) for the scalar product and 
its clasicial limit (4.11)-(4.2). Using this result, we were able to write a compact expression for the 
correlation function of three non-BPS operators in maximally supersymmetric Yang-Mills theory in 
the classical limit when the traces become large. Our starting point was Foda's determinant expression 
for the three-point structure constant [30]. The determinant formula of [30] was derived supposing 
that the three operators are deformed by a set of inhomogeneity parameters, whose values can be 
chosen at will. We computed the bosonized determinant expression for inhomogeneous problem and 
took, as in [30], the homogeneous limit at the very end in order to avoid spurious singularities. 

Another reason to treat the inhomogeneous problem is that this allows to extend, as argued in 
[56,57], the tree-level result to higher orders in A. Gromov and Vieira [56] showed that knowing 
the tree level solution for C^ 2 \ in presence of impurities, one can obtain the one-loop and the two- 
loop corrections by applying a special differential operator acting on the inhomogeneity parameters 
#(2) an( j #(3) Serban [57] proposed that this statement can be extended to all loops in the BDS 
model [58], i.e. when the dressing phase is not taken into account, and in the limit of large lengths 
Li,L2,L%. The analysis of [57] leads to the prescription that the higher loop corrections can be 
taken into account only by modifying the pseudo-momenta. For example, the three-loop result for 
the structure constant is obtained by changing the pseudo-momenta p u ,Pv and p w according to the 
three-loop Bethe ansatz equations [58]. 

According to [57], for finite value of the 't Hooft coupling A, the structure constant for three 
classical operators in the BDS model is given by (5.18) with a particular distribution of the Lj in- 
homogeneity parameters 9 in the interval [—2g,2g], where A = 16ir 2 g 2 . For this distribution the 
resolvents for the inhomogeneities associated with the three chains are given by 



This is the distribution of Li equal charges confined to the segment [— 2g, 2g] in absence of external potential. 
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where x is the 'Zhukovsky variable' defined by 



u = x + g 2 /x. 



(6.2) 



This has the same effect as changing the vacuum eigenvalues a(u) and d(u) of the transfer matrix, eq. 
(2.11), to 



In the strong coupling regime it is convenient to perform the change the variable (6.2) in the contour 
integrals in (5.18). We should mention here that the replacement (6.3) as a possible way to take into 
account (some of) the loop corrections was previously discussed in [10]. 

Assuming that this conjecture is correct, the expression (5.18) with the choice (6.1) for the inho- 
moheneities will give the all-loops result for the structure constant for the BDS model in which the 
dressing phase is neglected. On the other hand, the effect of the dressing phase [59], in the limit of 
large Li, L2, L3, is that the pseudo-momentum is modified as 



The fact that (5.18) depends only on the three quasi-momenta p u ,Pv,Pw and the resolvents G 9 (i) = 
Lid u log x invites one to consider the possibility that the result at finite 't Hooft coupling A is given in 
the classical limit again by (5.18), with the three pseudo-momenta modified according to (6.4). 

It is natural to expect that a systematic approach to the correlation functions of heavy operators 
in the full field strength multiplet of N = 4 SYM should be some extension of the algebraic curve 
method used in the spectral problem in [12-14]. 10 In the case of three long-trace operators, the struc- 
ture constant C123 is expected to be described by the ensemble of three algebraic curves, associated 
with the pseudo-momenta p u ,Pv,Pw In order to build a general algebraic curve formalism, one 
should learn how to compute more general inner products, at least in the classical limit. An inter- 
esting development in this direction was reported by Wheeler [61], who wrote determinant formulas 
analogous to (3.25) are obtained for the generalized model with SU(3) symmetry. 

Finally, it is certainly worth to try to adjust the expression (5.18) for the non-compact rank-one 
sectors of the full symmetry PSU(2, 2|4), such as SL(2,R) and SU(1, 1), where the integration 
contours should be placed along the real axis. 
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